UNITARY REPRESENTATIONS OF WAVELET GROUPS AND ENCODING OF 
ITERATED FUNCTION SYSTEMS IN SOLENOIDS. 

DORIN ERVIN DUTKAY, PALLE E.T. JORGENSEN, AND GABRIEL PICIOROAGA 

Abstract. For points in d real dimensions, we introduce a geometry for general digit sets. We introduce a 
positional number system where the basis for our representation is a fixed d by d matrix over Z. Our starting 
point is a given pair (A, 15) with the matrix A assumed expansive, and T) a chosen complete digit set, i.e., 
in bijective correspondence with the points in Z'^/A^Z'*. We give an explicit geometric representation and 
encoding with infinite words in letters from X>. We show that the attractor X{A^,'D) for an afBne Iterated 
Function System (IFS) based on (A,!)) is a set of fractions for our digital representation of points in R''. 
Moreover our positional "number representation" is spelled out in the form of an explicit IFS-encoding of 
a compact solenoid Sa associated with the pair (A, 2?). The intricate part fTheorem lG.lSI l is played by the 
cycles in Z'' for the initial (A,I')-IFS. Using these cycles we are able to write down formulas for the two 
maps which do the encoding as well as the decoding in our positional D-representation. 

We show how some wavelet representations can be realized on the solenoid, and on symbolic spaces. 
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1. Introduction 

Let A he a. d X d matrix over Z (the integers), and assume that its eigenvalues A satisfy |A| > 1. In 
particular A is assumed invertible. Let Z'^[A~^] be the associated discrete group obtained as an inductive 
limit 

^Z'^ ^Z'^ ^ ... 

i.e., 'L'^\A~'^\ = U^qA~'^Z'^, and let Sa denote the corresponding dual compact abelian group, where duality 
is in the sense of Pontryagin; Sa is a solenoid. 

Let be the transposed matrix, and let i : Z'^[A^^] R'' be the natural embedding (of groups), so i 
is a homomorphism of the respective additive groups. Now let i : —> Sa be the dual homomorphism. 
Hence i embeds M'' as a subgroup of the compact solenoid Sa ■ We will need this generality in an application 
to the analysis of wavelet multiresolutions. Here we use M.'^ = R'', i.e., M'' is its own Pontryagin dual. In 
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a special case of this construction, for d — 1, our embedding i corresponds to ergodic theoretic flows on 
compact spaces studied in a variety of contexts in dynamics. 

Note {{A^x) = aAiiix)), x E W^, where aA ■ Sa — * Sa is the endomorphism induced by A. 

Motivated by classical number theoretic problems in digital representations of fractions, we explore here 
symbolic representations of points in Sa- Our analysis uses an extension of George Mackey's semidirect prod- 
uct construction from representation theory (section 3) combined with a study of a family of combinatorial 
cycles (section 4.1). In section 6 we show that our encoding theorems (sections 4 and 5) apply to the construc- 
tion of wavelet multiresolutions, i.e., for generalized wavelet constructions where scaling in M.'^ corresponds 
to the matrix multiplication x Ax, x £ R'*; and the corresponding Z-action Z x M"* 9 {k, x) t-^ A'^x E M'*. 

Motivation from physics: The renormalization group. The idea of scale invariance is old. Its best 
know modern formulation in mathematics takes the form of iterated function systems (IFS), e.g., [HutSlj : 
and in physics it takes the form of a Renormalization group (RG) . But scaling arguments are commonplace 
in pure (e.g., wavelets) and applied mathematics; for example in attempts at explaining turbulence. 

The renormalization group makes its appearance in physics in different guises, often as a mathematical 
trick to get rid of the infinities, for example in quantum field theory, see e.g., |Fed87| IB0S91| . As a pure 
technique, it obtained maturity with for example the papers |Fed871 IFed941 IKW88[ IKW91| among others. 
The technique was developed also in quantum electrodynamics by R. Feynman, K. Wilson and others. The 
physicists devised theories of mass and charge renormalization. 

Old-style renormalization group (RG) techniques in physics have run into difficulties with non-renormalizability 
of gravity. Still they are used in various guises as tools in solid state physics, as they often get around di- 
vergence difficulties with the use of perturbation theory. 

As with iterated function systems (IFS), renormalization groups in physics attempt to describe infinite 
systems in terms of block variables, i.e.: some magnitudes which describe the average behavior of the various 
constituent blocks, often approximately true in practice, and good enough, to a first approximation. This 
is more or less amounts to finding long term behavior of a suitable RG transformation. When iterated 
many times, this RG transformation leads to a certain number of fixed points analogous to those seen in 
(non-contractive) IFSs. 

It was suggested in the 1970s by Don Knuth and others that there is an intriguing geometry behind 
computations in a positional number system. As is known since Euclid, representation of numbers in a fixed 
basis entails expansions in powers of the chosen base, say b. There is then a subset T> of the integers Z of 
cardinality |6| such that the corresponding "digital" expansion of real numbers is encoded by finite or infinite 
words in the "alphabet" V. In fact Knuth [Knu76 stresses that for a fixed b, there are many choices of digit 
sets V which yield a positional number system in a sense which is made precise. Similarly Knuth suggested 
the use of a matrix in place of b. In sections 5 and 6 below, we outline both the initial suggestion, the 
relevant literature, starting with [Odl78]; and we present our main results. The first sections of our paper 
address tools from representation theory central to our problem. 

Both the choice of the base b, and the set of digits allows for a great deal of freedom, even if we restrict 
numbers on the real line. Don Knuth suggested that this idea works in higher dimensions, i.e., in encoding 
points in M'^ this way. The case d ~ 2 oi course includes positional representations of complex numbers, and 
associated computer generated images in the plane. But by increasing the dimension further, this suggests 
using, for scaling basis, instead a fixed dhy d expansive matrix over Z, say A, in place of the base number 6, 
and a subset V of Z'^ of cardinality | det A\ for digits. This leads to the puzzling question: "In this geometric 
formulation, what are then the fractions as subsets of R'' ?" Early computer calculations by Knuth in 2D 
suggested that "fractions" take the form of "dragon-like" compact sets, "Twin Dragons" etc. Knuth's idea 
was taken up in later papers (by other authors, e.g., Lagarias and Wang) under the name affine Iterated 
Function Systems (IFSs) , and the set of fractions associated to a fixed pair {A, V) in d real dimensions have 
been made precise in the form of attractors for the IFS defined from (A,!?), see e.g., |Hut81j and B J99] . 

For points in Euclidian space, we introduce matrix scaling and digit sets. Our aim is to study the interplay 
between associated spectra and geometry. For a fixed matrix scaling and digit set we introduce a positional 
number system where the basis for our representation is a fixed d by d matrix A over Z. Specifically, a pair 
(j4, V) is given, the matrix A assumed expansive, and a finite set T) a chosen as a complete digit set, i.e., the 
points in V are in bijective correspondence with the finite group 11^ jAU^ . 

Such higher dimensional "number systems" allow more fiexibility than the classical one, introducing a 
computational device for the study of, for example, exotic tilings, wavelet sets and fractals. These are 
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geometric structures in studied recently in [BMM991 IChoOTl IJKS071 IShu03l IHutSll IBHS051 IFMM061 
lRud8 9]. We take advantage of a natural embedding of Mf^ in an associated solenoid, and we obtain an 
explicit solenoid-encoding of geometries in M'*, giving insight into notions of redundancy, and offering a 
computational tool. 

This expanded view also introduces novelties such as non-comniutativity into encoding. We give an 
explicit geometric representation and encoding for pairs pair (A,!)) (Theorem I4.9p . i.e., an encoding with 
specific infinite words in letters from V. 

Our positional "number representation" takes the form of an explicit IFS-encoding of points in a compact 
solenoid SA associated with the pair (AjT)). A crucial part (Theorem I6.15P is played by certain extreme 
cycles in the integer lattice Z"* for the initial {A, I?)-IFS. Using these cycles we write down a formula for the 
two maps which do the encoding as well as the decoding in our positional ^-representation. 

We will need basic tools from spectral theory (e.g., |Arv02j ) . but our aim is computational, still using 
operator algebraic tools, analogous to those used in the analysis of graphs and generalized multircsolutions, 
see e.g., [BMMSS, C hoQ7, JKS07|, IShu03| . Our use of iterated function systems follows conventions from 
|Hut81l IBH^ IFMMOfij . 

When an invertible matrix A and a finite subset T) in M.'^ are given, then we consider an associated finite 
set of affine mappings t^, indexed by points d in P, as follows, 

Td{x) := A-\x + d). 

Under suitable conditions on the pair {A, V) , see Definition 16. H repeated iterations of the combined system 
(jdjd^Tii then yields certain limit concepts. They take the precise form of either certain compact subsets in 
(attractors) or the form of limiting measures, so called equilibrium measures. In this formulation, the 
theory was made precise by J. Hutchinson in the paper |Hut81j . and this gave rise to what is now known 
as affine iterated function systems (IFSs). As is known, e.g.. |Hut81| . to make the limit notions precise, one 
introduces metrics on families of compact subsets in M'', or on families of probability measures. (In their 
primitive form these metrics generalize the known Hausdorff distance.) 

In harmonic analysis a subclass of the IFSs have been studied extensively by R. Strichartz and his co- 
authors, see |Str05[ IStr06| . The emphasis there is on discrete potential theory, while our present focus is on 
tiling and coding questions. However much of our motivation derives from harmonic analysis. 

It was also realized that wavelet algorithms can be put into this framework, and for fixed {A, 2?), one is led 
to ask for Haar wavelets, and to wavelet sets. In this paper we show that there is a representation theoretic 
framework for these constructions involving a certain discrete group that is studied in algebra under the 
name of the Baumslag-Solitar group. 

In addition to our identifying wavelet sets in a solenoid encoding (section 4), in section 6 we further 
show that the positional number representation for associated to a given pair (A, T)) takes an algorithmic 
form involving both "fractions" and "integer points". Here the fractions are represented by a Hutchinson 
attractor X{A^ ,1)) and the "(A, 2?)-integers" by a certain lattice F which makes X{A^ .V) tile by F 
translations. To compute this lattice F which makes an X{A"^ ^T>) tiling, we use a certain spectral duality 
(Lemma I6.25|) . Finally (section 6.2) the last mentioned duality is illustrated with specific planar examples, 
where both cycles and lattices are worked out. 



2. Iterated Function Systems 

The setting of this paper is a fixed d by d matrix A over the integers Z, satisfying a certain spectral 
condition, and its relation to the rank-c? lattice Z'* in W^. 

Traditional wavelet bases in L^(K'') are generated by a distinguished finite family of functions in L^(R'^) 
and the operations translation with the rank-d lattice Z'', and scaling with powers of A^ , j E Z. But there 
are similar wavelet constructions, super wavelets, in other Hilbert spaces which we explore here. See also 
}B.TMPn5l lDJ06bl ID,T07aj 

From the initial matrix A we form the discrete group 'l/\A^'^\ generated by powers of A\ j e Z, applied 
to Z''; and the compact dual solenoid group Sa- This solenoid is related to [0, 1)'' x where 51 is a compact 
infinite product of a fixed finite alphabet. But the two Sa and [0, 1)'* x fl are different, and their relationships 
are explored below. 
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First recall that matrix multiplication by A induces an automorphism a a in Sa- We are interested in 
periodic points for this action, and in a certain family of extreme orbits called cycles. As is well known, M'' 
is naturally embedded in Sa- 

Here we continue the study started in JDut06] of the support of the representations associated to super- 
wavelets. Starting with an embedding i of M'' into the solenoid, we have some periodic characters Xi in 'Sa, 
associated to the cycles. Then we show that the representation is supported on the union of where 
the multiplication here is just the multiplication in Sa- 

We combine Mackey's theory of induced representations with the analysis of iS^i-cycles. In this connection, 
we find a dynamical obstruction for embeddings as follows. Intuitively, one wants to encode the solenoid 
into a symbol space [0, l)'^ x 51. This can be a problem when one is dealing with matrices, as compared to 
the dyadic scaling in one dimension which is the traditional context for wavelet analysis. 

The reason for the obstruction is that the candidate [0, l)"* is not invariant under the inverse branches 
Td{x) = {A'^)~^{x + d), d ^ V where I? is a chosen finite set of vectors in R''. These maps serve as inverse 
branches to the action by on R'^/Z'^. The other candidate different from [0, 1)'' is the Hutchinson attractor 
X{A^,V) for Td, X(A^,V) contained in R'^. But X{A^ ,V) might not tile by Z'^. 

So one question is: Can one choose I? to be a complete set of representatives for the finite quotient group 
Z'^/A^Z'^ such that the attractor X{A^ ,V) of {ra} tiles R'^ by Z'^1 If not, then how should one choose A 
such that this is possible? Or perhaps, one must replace A with an iterate A^ of A for a suitable pi 

The reader may find our use of the pair of scaling matrix A, and its transposed A^ confusing. It is 
unavoidable and is dictated by our essential use of Fourier duality: If A is acting in d-space, then A"^ is 
acting in the dual vector variable (say frequency), see e.g.. Lemma 16.251 In deciding tiling properties for 
X = X{A^,D) our use of spectral theory is essential as there seem to be no direct way of attacking the 
tiling problem for X. 

While initially Knuth's analysis [Knu69| of what is now called affine Iterated Function Systems (IFS) was 
motivated by the desire to introduce geometry into algorithms for general digit sets in positional number 
systems, the idea of turning "digits" into geometry and tiling questions was followed up later by others, 
e.g., Odlyzko |Odl78| . Hutchinson |Hut81| . Bratteli-Jor gensen [BJ99] . and La garias-Wang [LW96a[ ILW96bl 
ILW96c[ ILW97[ ILWOO] . In [LW96aj the authors suggest that the tiling issues implied by the geometric 
positional "number" systems are directly connected with wavelet algorithms. In particular they pointed out 
that, for a fixed choice of A and V, the corresponding attractor X{A^ ,1)) as described above does not always 
tile R'' by translation vectors from the unit-grid lattice Z'^. When it does, we say that X{A^ ,1)) is a Haar 
wavelet. The term "Haar wavelet" is used because the corresponding indicator function is a scaling function 
(father function) for an ONB wavelet system in i^(R'^). The authors of |LW96a| showed that in 2D, every 
expansive d by d matrix over Z has at least one "digit" set V such that X{A'^ , V) is a Haar wavelet. It was 
later proved that in 5D, not every expansive 5 by 5 matrix over Z can be turned into a Haar wavelet; i.e., 
for such a matrix A, that there is no choice of T) for which X{A^ , V) is a Haar wavelet. 

The fact that there are exotic 5 by 5 expansive matrices A over Z, i.e., A in A^5(Z) for which no digit 
set V may be found such that X{A^ ,1)) makes a Z^-tiling of was worked out in the following papers 
|LW96cl ILW97[ IHLR021 IHL04j . By digit set we mean a subset T) in Z^, in bijective correspondence with 
Z^ /A^Z^ . Such exotic matrices A are said to not allow Haar wavelets. The question came up after Lagarias- 
Wang |LW97| showed that every expansive A in JVl2{Z) allows digit sets in Z'^ which make Z^ tile, i.e, they 
allow Haar wavelets. 

The aim of this paper is to revisit the geometry of sets X{A^ ^ V) in light of recent results on IFS involving 
dynamics and representation theory, see e.g., [Dut06[ [D J06al ID J06dl lDJ06c[ ID J07a| . 

3. Definitions and notations 

While standard wavelet bases built on wavelet filters and on a fixed expansive d by d matrix over Z, say 
A, refer to the Hilbert space LP'iW^), many naturally occurring wavelet filters [Dut06| suggest other Hilbert 
spaces than L'^{W^), in fact Hilbert spaces containing a copy of L'^{W^). This approach |Dut06j suggests the 
the name "super wavelets" , and naturally leads to representations of a Baumslag-Solitar group built on the 
matrix A, called wavelet representations. Starting with a fixed A, there is a compact solenoid Sa with the 
property that matrix multiplication by A on R'^/Z'^ induces an automorphism aA on Sa- 

In the past decade the literature on self-affine sets, encoding and digit representations for radix matrices 
has grown; in part because of applications to such areas as number theory, to dynamics, and to combinatorial 
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geometry. It is not possible here to give a complete list of these directions. Our present work has been 
motiva ted by the papers [XSOS] . |Cur06j . |GY06j . |HLR02j . [HL04] . |KLSW99j, ^Li06j . [LiOTj, [Saf98j . ^ZLZ06j . 
|LWOO| . 

Definition 3.1. Let Ahe a, d x d matrix with integer entries. We say that the matrix A is expansive if all 
its eigenvalues A satisfy |A| > 1. 

Let Ahe a d X d expansive matrix with integer entries. Let 

(3.1) Z'^[A-^] = {A-^k\j eN,k<EZ'^}. 
Note that Z'^[A~^] is the inductive limit of the group inclusions 

Z"^ ^ A-^Z"^ ^ A-^Z'^- ^ ... 

or equivalently 

z'^ 4 z"* 4 z'' 4 . . . 

On Z'^[A~^] we consider the discrete topology (even though H\A^^\ is a subgroup of 

On the group Z'^\A~'^\, the map aA{x) = Ax, x G Z'^[A~^] defines an automorphism of the group Z'^[A^^]. 

For the use of these groups in C*-algcbras, see e.g. |BJ99[ iBJKROlj . 

3.1. The group Ga Z'^[A^'^] xi^^j Z. The group Ga := Z'*[v4"i] x^^ Z is the semidirect product of 
gdj^-ij under the action of Z by the automorphisms a a- This means that 

(3.2) GA:^{{.l,b)\j eZ,beZ''[A-']}, {j,b) ■ {k,c) = {j + k,a\{c) + b), {j,k € Z,b,c E Z^iA-']). 

Proposition 3.2. The group Ga is generated by the elements u := (1,0) and tk — (0,fc), k E Z^[A^^]. 
Moreover 

(3.3) utkU-^^tAk, (fceZ) 

(3.4) t^-„fc := (0,^-"fc) = ij,A-''k) ^tA-nkU^, (n > 0, fc e Z^ i e Z). 



Remark 3.3. From Proposition 13.21 we infer that a unitary representation of the group Ga is completely 
determined by giving some unitary operators T^, k E Z'^ and U subject to the relation 

(3.5) UTkU~^=TAk, (fceZ). 

In |DJ07b| we use induced representations of Ga in the sense of Mackey in order to encode wavelet sets 
for a fixed expansive d by d matrix over Z. Note that Mackey's method was developed for continuous groups, 
where encoding is done with co-adjoint orbits. In contrast we show in the present paper that encoding in 
the solenoid is required for wavelet representations, i.e., representations of discrete versions of higher rank 
ax + b groups. 

3.2. The dual group of Z'^[A^^]: the solenoid Sa- The dual group of Z'' is T"*, where 

T'^ :={(zi,...,Zrf)||z,| = l,ze{l,...,d}}. 

For a; = (xi, . . . ,Xd) E M"* let 

For fc = (/ci, . . . , kd) E Z'^ and z = (e^^^^S . . . , e^"^'') E T"*, we use the notation 

^-Kikixi^ y2-KikfiX^ ^Trik-x ^ 

The duality between Z"* and T*^ is given by 

(3.6) (/^,z) = z^ {k eZ"^, zeT'^). 
For z = (e^'^"!, . . . e^'^"'') E T*^ we denote by 

(3 7) -yA ._ ^2TTi(A^)x _ ^^2TTiJ2j = iajiXj ^ ^ ^ ^ ^^2TTiJ2j = iajdXd-^_ 

Note that z'^'' = (z^)*^ for aU kEZ'^. 
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Definition 3.4. The dual group of Z'^[A ^] is the group Sa defined by 

Sa := {(z„)„gN I Zn e T'', = Zn, for aU n e Z}. 

The group Sa is called the solenoid of A. It is a compact Abelian group with the topology induced from the 
product topology on T^. 

The duality is given by 

(3.8) {A-^k , (z„)„eN) = {k , z,) = z], (j e N, A: G Il\ (z„)neN e 5^). 
The dual of the automorphism a a on Z'^[A^^], aA{x) = Ax is the shift 

(3.9) aA{zo,zi,...) = {z^,zo,zi,...), {{zq, zi, . . .) G Sa)- 

We denote by 0„ the projection maps 6'„ : Sa — > T^, 6'„(zo, zi, . . . ) = z„ for n e N. Note that 

(3.10) On+lO CTA = 071, i9n+lizo,Zi,...))^ = 9n{zo,Zi,...), (n € N, (zq , Zl ,...)) . 

4. EMBEDDINGS of R'' into the SOLENOID Sa 

We saw in Proposition 13.21 that there is a natural semidirect product discrete group Ga which carries 
a unitary wavelet representation. For wavelets in M.'^, this construction begins with a fixed expansive d 
by d matrix A, and the unitary representation will be acting in the Hilbert space L^(R'^). It is known 
that a certain redundancy ( |BDP05] . [HLOOj ) in wavelet constructions dictates unitary representations in 
Hilbert spaces larger than _L^(R'^), i.e., with L'^(R'^) embedded as an isomorphic copy in an ambient Hilbert 
space. By introducing a specific embedding of E** in an ambient solenoid Sa we are able to account for 
super representations fPefinition 14. 5p . The action of A induces an automorphism a a in Sa- By computing 
periodic points for a a we are able (Theorem 14. 9p to account for the super representations acting in an 
space defined from an induced measure on Sa- 

In sections 5 and 6 we will further study periodic points and cycles. It turns out that the notion of 
cycle is different when referring to the integer points and the fractions. Starting with a fixed radix pair 
(AjV), we will make "integer points" precise in terms of associated lattices (rank-d subgroups in E*^), and 
the "fractions" will take the form of compact subsets X in R'' defined by an {A, T>) self-similarity (Definition 

EH). 

We begin by showing how the space E'' can be seen as subspace of the solenoid Sa- 
Proposition 4.1. The inclusion i : Z'^[A^^] R'' has a dual i -.R"^ ^ Sa 

(4.1) l{x) = (e2-(-4")""-)„eN, {x G E-^). 

The map i is one-to-one, and onto the set of sequences {z„)neN G Sa with the property that lim„_>oo Zn — 1 
where 1 = (1, . . . , 1) is the neutral element o/T''. 
The map i satisfies the following relation 

(4.2) liA'^x) = aA{i{x)), (xeR"*). 

Proof- To see that i is one-to-one, we notice that if x, x' E R'^ and i{x) ~ i{x') then {A^)^'^x~ {A^)^'^x' £ Z"* 
for all n G N. Since A is expansive, the norm \\{A^)^'^''{x — x')\\ converges to as n ^ oo. Thus, x and x' 
must be equal. 

Since A is expansive, so is A^ , so {A^)~"'x converges to for all x G R''. Therefore e^^"^^ ^ converges 
to 1. Conversely, suppose {zn)n&fi is in Sa and z„ converges to 1. Then z„ = e^^"" for some Xn £ R'' and, 
for n large we can assume x„ is close to 0. Since we have z^_|_i = Zn, this implies that A^Xn-^-i = Xn modZ'', 
so A^Xn+i — Xn + I for some / G Z*^. But since both Xn+i and x„ are close to 0, this implies that, for 
some no G N, we have {A'^)xn+i = Xn for n> uq. Let x :— {A'^)"'''Xno- The previous argument shows that 
{A'^yxno = Xno-j modZ'^ for all j < rip, so i{x) — [e^^^^^ ^ "^)nen- Thus, i is onto. 

The other assertions follow from some direct computations, using the duality in (j3.8p . □ 

Now that we have the embedding of R'' into the solenoid, we can transport the wavelet representation on 
1/2 (M'^) to the solenoid Sa- 
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Definition 4.2. On L'^{R'^) we denote by Tk the translation operator {Tkf){x) = f{x - k), g £ L'^{W^), 
X S M'', k S Z"*, and by U the dilation operator {U f){x) = —j===f{A^^x). Their Fourier transform is 



(4.3) (ffc/i)(a;) =e2'^*'=-^/(x), {Uh)(x) = ^ \ del A\h{A^ x) , {h e L^{R'^),x eR'^,k e Z'^). 

The operators {U,Tk} (or {t/jTfe}) define a representation of the group Ga on i^(R'^). 

Definition 4.3. We denote by 5^(1) the set of sequences {zn)n&'i € Sa such that Um„^oo z„ = 1. On 
5yi(l) consider the measure fl defined by 

(4.4) / fdil= ^ /((z„)„eN)rfM(2), 



Where is the Haar measure on T''. 
On L^(5yi(l), /i) define the operators 

(4.5) ifkf)izo, zi, . . . ) = z^fizo, zi, . . . ), (zo, zi, . . . ) e 5^(1), fc G Z"*), 



(4.6) (C//)(zo,zi,...) = V|detA|/(aA(zo,zi,...)), ((zq, zi, . . . ) e 5^(1)). 

Our next theorem shows that this unitary resentation of the reduced solenoid 5yi(l) is a universal super 
representation of the wavelet group in that when A is given, then via an intertwining isometry W : L^(R'^) 
i^(iSA(l)) the standard A-wavelet-unitary representation acting on L^(K'^) is naturally included in the Sa{^) 
representation. The details of the symbolic encoding of this pair of representations depends on a choice of 
digit set T), and the structure of the associated (A, I?)-cycles (we will give these details in Section 6). 

Theorem 4.4. (i) The measure jl satisfies the following invariance property 



(4.7) 



(a) The operators T^, fc G Z'' and U are unitary and they satisfy the following relation 

(4.8) UniJ-^ = Tau: {k G ■L'') 

so {[/, Tfc} generate a representation of the group Ga on L"^ {S a{^) , fi) ■ 

(Hi) The map i is a measure preserving transformation between and 5^(1). 

(iv) The operator W : i^(R'^) L"^ {S a{^) , i^) defined by W/ = f o i~^ is an intertwining isometric 
isomorphism, 

(4.9) Wfk=fkW, (fcGZ'*), WU^UW. 

Proof, (i) The Haar measure on T'' satisfies the following strong invariance property: 

(4.10) / fdf^=T^f ^/(2/)dMW, {feL\f,)). 

Jfd det A J-fd ^ 

Using this, we have for / G L^(5^(l), fi): 

/ f°crAdjl=j /(z^,zo,zi,...)d^(z) = 

^ ' (z„)„ecSA(l),zo=z 

/" V V /(z^,zo,zi,...)d/x(z) = — ^— - /" V f{wo,Wl,...)d^l{z) = 

JTr-" ^ ^ _ \dctA\Jja .f-'.^. 



\delA\ J J. 



y*=z (z„)„e5A(l),2o=a («>„)„e5A(l),tuo=z 

f dfi. 

Sa{1) 



\ActA\ 

(ii) follows from (i) and some direct computations. 

(iii) First note that, by Proposition O (z„)„ G 5^(1) with zq = e^''" iff (z„)„ = i{y) and e^'^^^ = e^''*^, 
i.e., (z„)„ — i{x + k) for some k G Z'*. 
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Let / G £1(5^(1), A). Then 



/ f(i{x))dx^ I y2 f{i{x + k))dx = 
/ Y] f{zo,zi,...)dx = V f{zQ,zi,...)dfj.{x) = 



(^„)„e5A(l),^o=e2'"" (2„)„e5A(l),zo=z 

This proves that i is measure preserving. 

(iv) Since i is measure preserving, W is an isometric isomorphism. The intertwining relation (|4.9p follows 
by a direct computation that uses (|4.2p . □ 



4.1. Cycles. Next we will show how the "super-wavelet" representations from [BDP05| can be realized on 
the solenoid. 

Definition 4.5. An ordered set C := {Co, Cii ■ • ■ , Cp-i} IT'^ is called a cycle if = Q for j G {0, . . . ,p— 2} 

and C(f = Cp-i, where p > 1. The number p is called the period of the cycle if the points Q are distinct. 

Definition 4.6. Let C = {Co, Ci, ■ ■ • , Cp-i} be a cycle. Denote by 

(4.11) He ■■= L^(M'^) © • • • ® L^{W^) = L^(R'^ x Zp), 

. ^ 

p times 

where Zp = {0, . . . ,p — 1} is the cyclic group of order p. 
Define the operators on He- 

(4.12) TcAfo, • • • , fp-i) = (CoTfc/o, . . . , Cp~irfc/p-i), ((/o, . . . , fp-i) eUcke z") 

(4.13) Ucifo, . . . , /p-i) - (C//p-i, t//o, . . . , t//p_2), ((/o, . . . , /p-i) e Wc) 

where Tk and J7 are the operators on L^(R'') from Definition l4.2l We will denote by Tc^k and Uc the Fourier 
transform of these operators (the Fourier transform being applied on each component of Tic)- 

Then a simple calculations shows that: 

Proposition 4.7. The operators Tc\k, k E Z'^ and Uc are unitary and satisfy the following relation 

(4.14) UcTcmUc' = Tc,Ak, {k G Z'^) 
so they define a representation of the group Ga on Tic- 

Some examples and wavelet applications are also included in |Jor03| . 

Definition 4.8. Let C {Co, • • • , Cp-i} be a cycle, Q = e'^"''"' , for some Xj G M'', (j e {0, . . . ,p - 1}). We 
will use the notation On '■= dnmodp, Cn Cnmodp for all n €z Z. We denote by 

Xc (Co, • • ■ , Cp-i, Co, • • • , Cp-i, ■ • • ) e Sa- 

Let 

p-i 

(4.15) 5a(C) := |Ja7(xc)*(M''). 
Let ic : M.'^ X Zp ^ Sa{C) 

(4.16) lc{x,j) = o-j,\xcW) = (e2-^«^")""-+«"+^))„eN, {x G j G Zp). 

Define the measure fi on Sa{C) by an equation similar to (j4.4p . (the only difference here is the support 
<S^(C) instead of<SA(l)) 

/ fdjl^j ^ f{{zn)nen)dfi{z). 

JSa{C) JT-i (3„)„e5^(c),eo((^„)„eP,)=z 

Define the operators Tk, fc G Z'' and U on L^(5a(C), /i) by the same formulas as in (|4.5p and (|4.6p . 
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Theorem 4.9. (i) The point xc is periodic for ua, o'^(xc) = XC OLf^d ua permutes cyclically the sets 

(a) The set Sa{C) consists of exactly the points {zn)neN with the property that the distance from Zn to C 
converges to 0. 

(Hi) Let UA.p : Z'' X Zp ^ Z'^ X Zp 

(4.17) aAAx,j)^iA'^x,j~l), {xeR'^.j eZp). 

The map ic is a bijective measure preserving transformation that satifies the relation 

(4.18) ic °OiA,p = cTAoic- 
fiv) The operator Wc ■ He ^ i^(K'^ x Zp) L"^ {S a{C) , jl) , Wcf = / o , is an intertwining isometric 

isomorphism: 

(4.19) Wcfck^nWc, (fceZ-*), WcUc^UWc- 
Proof, (i) is trivial. 

(ii) If {zn)n is in SAiC) then for some j, we have that {cr^'^ {xc))~^{zn)n is in i(R'^) so C,~^Znp converges 
to 0. Therefore Znp converges to Cj so Znp-i = z^^ converges to Q^i for all I £ Zp. 

Conversely, suppose {zn)n € Sa and dist(z„,C') converges to 0. We claim that z„p converges to one of 
the points of the cycle Q. 

Pick an e > small enough such that for all i € 7jp, dist(z„,^i) < e implies dist(z„,Ci') > e for i ^ i' . 
There exists a (5 > 0, J < e such that for all i G Zp, if dist(z, Ci) < <^ then dist(z'^'', Ci) < e. 

There exists an rie such that if n > then dist(zn,C) < ^ < e. Then for some i £ Zp we have 
dist(z„^,Ci) < £• Also for some i' € Zp we have dist(2;„^+p, Ci') < <5. This imphes that dist(z„^, ^i') < e so 
i' = i. By induction, we obtain that dist(z„^+fcp, C.i) < S < e. And since dist(z„, C) converges to this shows 
that dist{zn^+kp, Ci) converges to 0. Applying the map z^ several times we conclude that Znp converges to 
one of the elemenst of the cycle, Q. 

Then consider (u;„)„ := (JAixc)^^{zn)n & Sa- Clearly Wn converges to 1. By Proposition 14. 1[ there 
exists an X gM."^ such that {Wn)n — i{x). Thus (z„)„ = a''j^{xc)i{x) j and this proves (ii). 

(iii) Since i is bijective (Proposition 14. ip . clearly ic is also bijective. To check that ic is a measure 
preserving transformation, since i is measure preserving by Proposition 14. 1[ it is enough to check that 
multiplication by cr^"'(xc) leaves the measure jl invariant, i.e., for a function / defined on (t^-'(R''), 



(4.20) / fdfL= f{(JA'{xc){zn)n)dfL{znU (jeZp). 

It is enough to check this for j = 0. Using the translation invariance of the Haar measure /i on T*^, we have: 
/ fdfi^ I V /((z„)„)d^(z) / V f{xc{wn)n)d^i{z) = 



Equation (|4.18p follows by a direct computation. 

(iv) Follows from (|4l8)) . □ 



X! f{Xc{Wn)n)d^l= f{Xc{Wn)n)dfl{Wn)n 



5. Encoding of integer points 

The idea of using matrices and geometry in creating a positional number system for points in was 
initiated by Don Knuth, see especially |Knu76| . vol 2, chapter 4 (Arithmetic), and 4.1 which introduces this 
geometric and algorithmic approach to positional number system. In fact, the Twin-Dragon appears on page 
206 (in V 2 of |Knu76j ) . In our present discussion, with a fixed expansive matrix playing the role of the 
basis-number (or the radix) in our radix representations, the natural question arises: "What is the role of 
the integer lattice Z'^ relative to our radix system?" This section gives a preliminary answer to the question, 
and the next section is a complete analysis involving cycles. 
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As before, we begin with a choice of expansive matrix A (i.e., A is a fixed d hy d matrix over Z), and 
we choose a subset T) in Z'' (for digits), points in V in bijective correspondence with Z"^ / Z'^ . As it turns 
out, "the integers" relative to the (^,I?)-radix typicaUy wiU not have finite radix (or Laurent) expansions 
in positive powers of . The reason for this is the presence of certain non-trivial cycles C in Z'* leading to 
infinite repetitions, which we will take up systematically in the next section. 

In this section we will show that when the pair {A, V) is fixed, there is an encoding mapping which 
records the finite words in the alphabet D which will correspond to the cycles in Z'' that are associated with 
our particular choice of (A,!)). However, once the cycles in Z'* are identified, there is a much more detailed 
encoding directly for Z'' which will be done in detail in Theorem 16. 151 

So our present encoding mapping (j> : ^ I?^, depending on the pair (A, I?), is an introduction to our 
analysis of a refined solenoid encoding and of all cycles in the next section. Still our starting point is a fixed 
radix-pair (A, P). The fact that the encoding mapping (f) : ^ T)^ is injective is a consequence of the 
expansive property. Corresponding to {A, V) there is a finite set of finite words F = F{A, V) in letters from 
v. These words label the Z'^-cycles, and the encoding mapping </) : Z'^ ^ (infinite Cartesian product) 
maps onto the set of infinite words which terminate in an infinite repetition of one of the words from F. 

Let d be given. Let B he a, d x d matrix over Z, and assume 

(5.1) Pi B'^Z'^ = {0}. 

fc>i 

Note that this holds if B is assumed expansive. 

Let P C Z*^ be a complete set of representatives for Z'^/BZ''. Assume £ V. (This assumption is for 
convenience and can be easily removed mutatis mutandis.) 

Definition 5.1. Define a {B,!)) encoding of Z'* 

<j>{x) ■.= dodid2...,<l>:Z'' -^V^^ 

as follows: when x £ Z'^ is given, there is a unique pair xi £ Z"^ and do £ T) such that 

(5.2) x = do + Bxi. 

By the same argument, now determine di,d2, - ■ ■ £ T), and X2, X3, • • • £ Z''' recursively such that 

(5.3) Xk^dk+Bxk+i- 
Definition 5.2. Set 

00 

n:=v^ = Y[v. 

Elements lu £ CI are called infinite words in the alphabet V. If w is a finite word, we denote by v the infinite 
repetition of this word vvv .... If there are finite words v and w such that lu = (vw), we say that co ends in 
a cycle. 

Proposition 5.3. (i) The encoding mapping : Z'' — s- is well defined, 
(ii) (p is one-to-one. 

(Hi) (j) maps onto a subset of of all infinite words that end in cycles. 



Proof. Part (i) is immediate from (15. 2p . (15. 3|) . 

(ii) Suppose x,y £ Z"^ and (f>{x) — (f>{y). Then an application of (j5.3|l shows that x — y £ DkyiB'^Z'^, and 
we conclude that x = y hy an application of (|5.ip . 

(iii) Follows immediately from Thcorcm l6.15l □ 

Remark 5.4. (i) For x £ Z'', the encoding (f)(x) = vw ~ dodid2 . . . with di £ V, and v,w finite words is 
unique; but the formal sum 

(5.4) do + Bdi + B'^d2 -(-... 

is not convergent unless u; = = 000 . . . infinite repetition. In that case there exists m e N such that 
V = do . . . dm-i and 

X = do + Sdi + • • • + 
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(ii) Suppose ?; = and (f>{x) = w, with w = IqIi . . . Ip-i- Then —x has the following infinite, convergent, 
periodic fractional expansion 

oo p— 1 

See Proposition l6.14l and Theorem 16 . 1 51 for the proof. 

Example 5.5. The following simple example in ID illustrates the cases (i) and (ii) above. 
Let d = 1 , B = 2, and V = {0, 3}. Let </> : Z ^ be the encoding. 
We have 0(11) = 300330, i.e., v = 3003 and w = 30. 

And 4>{18) — 0330, i.e., v = 033 and w — 0, corresponding to the finite representation 

18 = + 3 -2 + 3 -22. 

Finally 0(— 2) = 03, i.e., v = 0, and w — 03. Hence by (ii) we get the following infinite fractional dyadic 
representation 

2 = 3 • 2^1 + • 2^2 ^ 3 • 2"'"^ + • 2"^ + 3 • 2"'^ + . . . 

Proposition 16.141 shows that the cycles are obtained by intersecting Z with the set —X{B,'D), where 
X{B,'D) is the attractor of the maps to{x) = x/2, T3{x) = {x + 3)/2. In our example X{B,'D) = [0,3]. 
There are two cycles of length one {0} and {—3} and one cycle of length two {—1, —2}. 

The encoding mapping records the cycles as follows: the one-cycles 0(0) — 0, 0(— 3) = 3. The two-cycle: 
0(-l) = 30, 0(-2) =03. 

The next section takes up the encodings in general. 

6. Encodings of the solenoid 

In this section we return to the geometry of general digit sets in positional number systems, turning 
"digits" into geometry and tilings. The starting point is a given pair (A, V) with A expansive over Z, and 
I? a complete digit set. With the aid of the solenoid we give an explicit encoding. Specifically, we show that 
the attractor X{A^,'D) for the corresponding afhne Iterated Function System (IFS) is a set of fractions for 
an X')-digital representation of points in E**. Moreover our positional "number representation" is spelled 
out in the form of an explicit IFS encoding of the compact solenoid Sa associated with the pair {A,!)). 
The intricate part (Theorem 16. ISp is played by the cycles in Z'^ for the initial {A,'D)-IFS. Using the cycles 
we are able to write down formulas for the two maps which do the encoding as well as the decoding in our 
positional ^-representation. 

Take a point {zn)n&'N in the solenoid Sa- Then zq G T*^, zf — zq, z^ = zi, and so on. 

Since zq is in T'^ it can represented by e^'^""^ where xq G M''. Note that one has several choices for xq, 
any of its integer translates xq + k, k E Z"*, will do. 

Then z^ — zq, so zi is a "root" of zq- There are |detA| choices: if zq — e^'^'^" then zi must be one 
of the points e^^'^'^^) ^{xo+d)^ d E V, where V is some complete set of representatives for l/' /A^U^. Say 

— g2iri(A^)^^(a;o+(io) „ g2iria;i 

At the next step Z2 is a root of zi so Z2 = e^'^*'^"^^) ^{xx+dx) some c?i e V. By induction, we get a 
sequence di, ^2, ■ ■ • in P. 

Thus, picking a point (z„)„gN in Sa amounts to choosing an xq e and an infinite word d^di • • • e O := 
I?'^. Thus we say that (2;„)„gN can be encoded as 

{Zn)nm ^ {X{),dodi .. . ). 

Now note that changing the choice of xq (to say xo + k), affects the entire sequence do, di, ■ • ■ • We want 
to make this choice unique in some sense. For this we need to find a subset F of M'', such that for each 
z € T'^, there is a unique x G F such that z = g^iria;^ other words, F must tile M'^ by Z'^-translations. 

Of course a first choice of this set F would be [0, 1)''. While this works in dimension d = 1, it may be 
inappropriate for higher dimensions. The problem is that we would need also zi to come from e^'^'^^ with 
xi e F, and this would mean that xi = {A'^)'~^{xo + do) is in F. Thus, our set F must have the following 
property 

\J{A^r\F + d)cF 

deT> 
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But since (zi, Z2, ■ ■ ■) is also an element oi Sa and zi can be any point in T'*, it follows that we must actually 
have 

(6.1) [jiA'^)-\F + d)^F. 

Of course, when we are interested only in measure theoretic notions, we can allow the equalities to hold only 
almost everywhere. 

When F is compact, this equation identifies F as the attractor of an afhne iterated function system. 

Definition 6.1. Let I? be a complete set of representatives for Z'^/A'^Z'^. For every d e I?, we denote by Td 
the map on R'^ defined by 

(6.2) Td{x)^iA^)-\x + d), (xeK'^). 

With this notation, if F is compact, equation (|6.ip says that F is the attractor of the afSne iterated 
function system {Td)de'D- This identifies F as 

oo 

(6.3) F = X{A^,V) {Y,{A^)-^d, \ d, G V}. 

i=i 

To find an encoding of the solenoid Sa means to find a subset F of Mf^ and a complete set of representatives 

V of Z'^/A^Z'^ such that, if x € F then r^x e F for aU d € V, and the decoding map d : F x V^* ^ Sa 
defined by 

FXV^3 {Xa,dadi ...)^ (g2,r»xo^ g27r»r,„xo^ g2^»r,^r,„:.o^ . . . ) G 5^ 

is a bijection. 

Thus, to find this encoding of Sa we need a subset F of M.'^ and a complete set of representatives V that 
satisfy (|6.ip and such that F tiles M'' by integer translations, i.e., 

(6.4) (J (i^ + fc) = K'' Lebesgue-a.e., and (F + fc) n (F + k') = Lebesgue-a.e. 

So the problem of encoding the solenoid into the space F x is equivalent to the following: 
Question. Given an expansive dxd integer matrix A, is it possible to find a complete set of representatives 

V of Z'^/A^Z'^ such that the attractor X{A'^,V) of the iterated function system lTd)dev tiles R'^ by Z'^7 
As explained in the introduction, while this is known to be true for dimension d = 1 or d = 2, there are 

some 5 by 5 matrices for which such a V does not exist. 

Definition 6.2. (a) We say that a subset F C K.'' tiles K*^ by a lattice F iff the following two properties 
hold: 

R'^= U(F + 7) 

(F + 7)n(F + 7') = 0, 7,7'er,7^7'. 

We say that F tiles R'' up to measure zero by the lattice F if these two properties hold up to Lebesgue 
measure zero. 

(b) By a lattice we mean a rank-d subgroup of R"*. We shall be interested in sublattices F C Z''. For a 
fixed F C Z'' we say that F is of index k if the order of the quotient Z'^/T is k. 

Lemma 6.3. Suppose a relatively compact subset F C R'' tiles by some lattice F C Z'^. Then the lattice F 
is of index k if and only if the mapping F 9 x i— > g^'^*^ G T'' is k-to-1 (up to measure zero). 

Proof. It follows from the definition that F tiles by F iff the restriction to F of the quotient mapping 
R"^ R'^/F is bijective up to measure zero. Hence the assertion that the given map is fc-to-1 is equivalent 
to the natural mapping R^'/F R''/Z'^ being a fc-fold cover; but this is so by the induced isomorphism 

(RVr)/(R''/Z'^) = Z'^/F. □ 

Proposition 6.4. Let fl :— . Define the map d : R''- x D, ^ Sa by 

(6.5) 0(x, iuouji . . . ) = (e^"", e2---o-, e2"'"-i"-o-, . . . ), (x G R^ luquji ■■■en). 
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(i) For each x £ W^, k e Z'', and cu G ft there is a unique ui' — uj'{x,k,uj) G such that (){x,uj) = 
d{x + k,uj'). Moreover, if x,x' S M.'^ and uj,uj' G fl, such that d{x,uj) — (3{x',uj'), then x' — x + k 
for some fc G Z"* and uj' = uj'{x, k, oj). 

(ii) Let F be a subset o/K'^. The restriction of the map'it to Fxil is injective if and only if Fn{F+k) = 
for all keZ'^, k^O. 

(iii) The restriction of Xl to F x is onto if and only if 

y (F + k) ^w^. 

(iv) The restriction of d to F x ft is bijective if and only if F tiles by Z'^ -translations. 

(v) Define the map p : R'^ x Q R'^ x Q 

(6.6) p{x,ujouJi . . .) = {T^gx,ujiuj2 ■ ■ {x eR'^^ujQUJi ■ ■ ■ e il). 
Then 

(6.7) do p = cj^^ od. 

Proof {i)Wewa,ntT^gX = T^'^{x+k). So {A^)~'^uJo = {A^)~'^{k+uj'Q), i.e., uj'q = oJo-kuiodA'^Z'^. Since I? is 
a complete set of representatives for Z'^/A'^Z'^, there is a unique ujq G T) such that this is satisfied. Proceeding 
by induction we see that u}'2,uj'^, . . . can be uniquely constructed such that e^'^^'^'^i-'-'^'^o^ = g27rir^/^...r^/^(3;+fe) 
for ah n e N. 

If ii{x,u!) ~ d{x',uj') then e^^" — e^^" so x' — x + k for some k E Z"^. The rest follows from the 
uniqueness of uj'{x, k, iv). 

(ii) Suppose restricted to -F x f2 is injective. Then, by (i), we cannot have x and x + k in F for some 
^ 0. Conversely, if is not injective on this set, then 5(x,(jj) = ^{x'jUj') for some x,x' E F and uj,lu' G fl. 

Using (i) again we get that x' = x + I for some I G Z"^ so F O {F + I) ^ (d. 

(iii) Suppose the restriction of to F x is onto. Then for all y G M'^, there is a; € F and ui G such 
that d{x, uj) — i{y). Then e^'^" = e^'^''' so y — x + k for some k G Z'', and therefore y G F + k. This shows 
that U{F + k) =M'*. 

Conversely, take (2n)neN S Sa- There exist a;„ G M'* such that z„ — e^^"" for all n. By hypothesis we can 
take xq G F. Since — zq we have that, A^xi = xq so for some uiq G T>, r^gXo = xi. Then, by induction 
we can construct lj„ e I? such that t^^^ . . .t^^xo = Xn+i- This proves that X){x,lo) = (e^'^"")„ = (zn)n- 

(iv) follows directly from (ii) and (iii). 

(v) requires nothing more that a simple computation. □ 

Proposition 6.5. Suppose F is a subset ofM."^ that tiles R'^ by a sublattice T ofZ'^ with \Z'^ /T\ — N . Then 
the restriction of the map d to F x is N-to-1. 

Proof. Since UkeZ'^i^ + k) D U^^r{F + 7) = R'^, it follows from Proposition 16.4( 111) that the map is onto. 

We claim that for each x G F there are exactly N points k G Z'^ such that x + k G F. Indeed, let di, . . . dN 
be a complete list of representatives for Z'^/F. Then for each i G {1, . . . , N} there is a unique G T such 
that X + di G F + "fi. Then we can not have di — ji = di' — 7^' for i ^ i' (that would imply di = di' modP) 
so the points di — 7^ are distinct and x + (di — ji) G F . 

Now we can use Proposition 16.4( 1) to see that d restricted to F x is A^-to-l. □ 

Consider now the compact attractor X{A^,T>) of the iterated function system (rd)de'D given in (|6.3p . It 
is known (see |LW96c[ ILW97| ) that X{A^,V) always tiles R'^ (up to measure zero) by some sublattice P of 
Z''. 

The connection between lattice tilings and spectral theory was studied systematically in |Fug74| and 
[Ped96| . and we will introduce spectrum in subsection 6.2 below. From the choice of digit set "D for a fixed 
matrix A, we conclude that X{A^ ,T>) has non-empty interior. In fact the d-dimensional Lebesgue measure 
of X{A^, V) must be an integer. It is 1 if and only if X{A^ , V) tiles M'* by the "unit-lattice" Z'^. By spectral 
theory we are referring to the Hilbert space L?'{X{A^ ,T>)). 

Definition 6.6. We say that {A,T>) satisfy the tiling condition if X{A^ ,T>) tiles R'^ (up to measure zero) 
by the lattice Z'^. 
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In subsection 6.2 below we give examples for d = 2 of pairs {A, V) which do not satisfy the tiling condition. 
Nonetheless, even if some particular pair (A, T)) in the plane does not satisfy the tiling condition, it will be 
possible to change the digit set T) into a different one V , while keeping the matrix A fixed, such that the 
modified pair {A, V) will satisfy the tiling condition. But by going to higher dimensions (d = 5) as we noted 
there are matrices A for which no V may be chosen making {A, V) satisfy the tiling condition. 

Lemma 6.7. For all d, d! T) , d ^ d' the intersection Trf(X(^"^, D)) n t^' (X{A'^ , T))) has Lebesgue measure 
zero. 

Proof. We have the following relations, with /i the Lebesgue measure: 

X{A^,V)= y Td{X{A^,V)), 



^,{Ta{X{A^,V))) ^ 



As a result we get 



H{X{A^,V)) = ^ ^ j i{X{A^,V)) - fi{ combined overlap sets). 
|detA| 

Therefore the combined overlap sets must have measure zero. □ 

Proposition 6.8. Suppose {A^T)) satisfy the tiling condition. Then the function p defined in (16. 6p maps 
X{A^ ,1?) X n onto itself and the restriction of p to X{A^ ,1)) x ft is injective a.e. in the sense that the set 
of points X G X{A"^ ^T)) with the property that there exist x' G X{A'^ ,1)), lu,uj' E p{x,uj) ~ p{x',uj'), has 
Lebesgue measure zero. The inverse p^^ of this restriction is defined by 

(6.8) p^^{x,LiJi)UJl . . .) = {A^ X — L0x,UJx^f3^1 ■ ■ ■), (x G I?), Wq'^I ■ • ■ S ^^): 

where lOx is the unique element ofT> with the property x G t^^{X{A"^ ,T>)) . 
Proof. Since 

(6.9) X{A^, V) = \Jdevrd{X{A^ , V)), 

it follows that p maps X{A^,V) onto itself. 

Suppose now p{x,uj) = p{x',uj') for x,x' G X{A'^,T>), uj,uj' G ft, and {x,uj) ^ {x',uj'). So either x ^ x' 
or a; 7^ uj' . When uj ^ lo' , since p{x,uj) — p{x' ,uj') it follows that ijJiuj2 ■ ■ ■ = uj[lu2 ... so cjq 7^ '^o- Also 
T^oX — T^'^x' . But T^„{X{A'^ ,T>))nT^'^{X{A'^ ,T>)) has measure zero (see Lemma l6.7p . so a; must be in a set 
of measure zero, li lo = lu' then loq — lo'q so Tj^qX — t^i^x' implies x = x' . This proves the injectivity of p. 

Since Td{X{A^ ,1))) are mutually disjoint, the element lo^ of T> is well defined. Then it is easy to check 
that p^^ is indeed the inverse of the restriction of p. □ 

6.1. Encodings of cyclic paths. Let I? be a complete set of representatives of jX^U^ . And let := Tf^ . 

Consider now a cycle C :— {Co, Ci, ■ • ■ , Cp-i} and suppose Q — e^^*^^ for some dj G K*^. If in addition, A 
and T> satisfy the tiling condition we can pick 6j in X{A'^ , T>). Then, since = zq, there is a G I? such that 
6i — ti^Oq. Continuing this process, we can find Iq, li, . . . , Ip^i such that Ti.Oj = 0j+i for j G {0, . . . ,p — 2} 
and Ti^_^9p^i = Oq. Thus the point xc = (Co, • ■ • , Cp-i, Co, • ■ • , Cp-i, Co, ■ • ■ ) from Sa can be encoded by 
(^0, ^0 • • ■ Ip-ila ■ ■ ■ Ip-i, la ■ ■ i-e., by an infinite repetition of the finite word Iq . . . Ip^i. 

Definition 6.9. A finite set in M'', C — {Oq, . . . , 9p^i} is called a cycle if there exist Iq, . . . , Ip^i G T> such 
that TigOo = 6*1, Til 6*1 = 02, . . . ,Ti^_20p-2 = dp-i and Ti^_^9p-i = Oq. Thus 6*0 is the fixed point of ti^_^ ■ ■ .n^, 
9i is the fixed point of ti^ti^_^ ■■ - Ti^, . . . , 9p-i is the fixed point of ti^_2 . . . n^n^^^. 

The points 9o, . . . , 9p^i are called cyclic points. We say that is the cyclic point associated to Iq . . . Ip-i, 
and we say that C — {0o, . . . , Op-i} is the cycle associated to Iq . . . Ip-i. 

Let C := {e^'^*^", . . . , e^'^^^p-i} be a cycle associated to Iq . . . Ip^i. Take now a point (z„)„eN in Sa{C). 
We want to see how the encodings of points in Sa (C) look like. 
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By Theorem 14. 9[ the point (z„)„gN is in one of the sets cr^"' (xc)*(R'^)- Suppose Zq = e^^" for some 
X € M'^. Then {zn)neti = o-a'' ixc)'i{y) for some y £ R'^, and fooking at the position, x = dj + y so 
y — X + k — 9j for some fc £ Z'' . Thus 

izn)neN = o-A^{xc)iix + k - 0j) = ic{x + k- 0jj). 
On the other hand, according to the previous discussion, (2;„)„gN is equal to (e^'^*^, e^'^^'^'^o^, g^-n-iT^-^r^gX ^ ^ ^ 
for some infinite word luquji .... Thus we must have some precise correspondence between the pair (fc, j) € 
Z'' X Zp and the infinite word ojquji ■ ■ ■ G ft. Since z„ is approaching the cycle C as n — > cx) one might expect 
that the infinite word ujqlui . . . ends in a repetition of the finite word ■ • ■ Ip-i that generates the cycle C. 
While this is often true, there might be some other cycles C that are congruent mod Z'' to C, that will affect 
this encoding uj. In any case, ujqUJi . . . that corresponds to (fc, j) will be eventually periodic, and it will end 
in an infinite repetition of a finite word that corresponds to such a cycle C . 

Definition 6.10. We denote by ■ • ■ Ip-i the infinite word in fl obtained by the infinite repetition of the 
word lo ■ ■ ■ Ip-i- Let 

ilc {wo • • ■ cj„ ^o . . . Ip^i \ ujo, . . . ,u!n e V, n e N}, 
i.e., the set of infinite words that end in an infinite repetition of the word lo . . . Ip^i. 

There are some cycles which have cycle points that differ by integers. Such cycles would make our encoding 
ambiguous, so we avoid this situation. 

Example 6.11. Let d = 1, A = 2 and V = {0, 3}. Then tqX = x/2, T3 = (a; + 3)/2. Then it is easy to check 
that the attractor X{A^ ,!)) is [0,3]. The set {1,2} is a cycle that corresponds to 30, and its points differ 
by an integer. 

Definition 6.12. We say that the cycle C = {9o, . . . , dp-i} is simple if 9j ^ Ojr modZ'* for j ^ j' . 

Following |BJ99] . for a simple cycle C = {6*0, . . . , Op-i}, we define an automorphism TZc on the set Z*^ — C. 
Note that since the cycle is simple, the sets Z'* — 6j are mutually disjoint. The map TZc is an extension of 
the division with remainder. Here we "divide" by A'^ . For each point in a — 6j G Z'' — 6j, there is a unique 
"quotient" b — 0j+i in Z'' — 0j+i and a unique "remainder" do G D such that 

a-0j ^ A'^{b - 0j+i) + do. 

Then TZc {a — Sj) is defined as the quotient TZc {a ^ Sj) — b — 0j+i- 

We used here the fact that A'^9j = 0j-i modZ'', because Ti._^0j-i — 0j, for all j e Z. Recall also that 
we use the notation 0j := 0jmodp for j G Z. 

Definition 6.13. Let C = {0o, ■ ■ ■ , 0p-i} be a simple cycle. On Z"^ — C = Uj=o ^ ^j) define the map 
TZc as follows: for each a e Z'' and j £ {0, . . . ,p — 1} there exist a unique b £ and do £ T> such that 

(6.10) a-0^ = A^ib - 0j+i) + do. We define 7^c(a - 9j) := b - 0j+i. 

Therefore TZc (a ~ ) is defined by 

(a - 0j) - A^TZc{a - 0j) £ V. 

Also, ~TZc{a~0j) = Td„ (— (a — 6'j)), where do is the unique element of V such that T^g {—{a~0j)) £ Z'' — 6'j+i. 

The encoding of (fc, j) e Z'' x Zp is obtained by a generalized Euclidean algorithm: take k — 0j in Z''' — 0j, 
then "divide" by A^ and keep the remainder: k—0j — A^TZc{k~0j)+ujo. Then take the quotient TZc{k-'0j), 
divide by and keep the remainder uji , and so on to infinity. The infinite sequence of remainders will give 
us to. 

But first, we need some properties of the map TZc. 

Proposition 6.14. Let C — {0o, . . . ,0p-i} be a simple cycle. Let X{A'^,T>) be the attractor of the iterated 
function system {Td)de'D. 

(i) A point t £ C — 1^ is a cycle point for the iterated function system {Td)de'D if o,i^d only if there is 
some n > 1 such that TZ'^{—t) — —t, i.e., ~t is a periodic point for TZc. Moreover if t is associated 
to Too . . . TOg_i then q is a multiple of p and 

m„ = 7^;^(-^)-^^7^J+l(-^), (neN). 
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(ii) For every t & C — 21^ there exists a I > such that is periodic for TZc, i-e., every point in 
1/^ — C is eventually periodic for TZc- Moreover —TZ'(j{—t) is in (C — Z'') n X{A^ ,1)). 

(iii) The intersection (C — Z'') H X{A'^ ,1)) consists exactly of negative the periodic points for TZc- 

Proof (i) lfto = te 6j - Z'* is a cyclic point for {Td)dev, then r„(,io = h, t„iJi i2, • ■ • , Tm,_iiq-i = to 
for some mo, . . . , mg_i e T> and some ii, . . . , iq-i G M"*. Then = ^"^to — so G 9j^i — Z"*, 

(because A^6'j = By induction e dj+i-q — Z'' for ^ € {g — 1, g — 2, . . . , 0}. 

Since to G ^ and also to G ^j-g — ^'^j as the cycle is simple, it follows that q must be a multiple of p. 

We have Tmoto — ti so —to = ^"^("^i) + "^o- Also — ti G Z'' — ^j+i- This means that TZc{—to) = — ti 
and mo = —to — ^"^7?.c(~to). By induction TZc{—tn) = — t„+i and m„ = — 1„ — A^TZci—tn)- This proves 
one direction. 

For the converse, if 7?.p(— to) = —to, for some to G 9j—Z'^, then for each n there is some to„ G I? such that: 
7^'^(-to) = ^^7^;^>+^(-to)+m„. Thus the sequence {to„} has period g, and t™„ (-7^J^(-to)) = -7^p+^(-to), 
which proves that to is in the cycle {—(—to), —TZc{—to), . . . , — 7^^^^(— to)}. 

(ii) Since A is expansive there is a norm on such that for some < c < 1, || (A-'")^^^!! < c||a;|| for all 
X G R'^- Then if i? > £i£«LM ^ 

t4B{0,R)) C B(0,i?). 

Indeed ||Tdx|i < c(||a;|| + ||d||) < c(i?+ |id||) < R for all x G B{0,R) and all d G P. 

Take now a ~ 0j G Z'^. Take some R > max{||a - 6*^ 1 1 , -^^^^i^^^M } . Then note that 7^c(a - 0^) = 
-rd(-(a-6'j)) forsomedG V. Therefore 7^c maps 5(0, i?) nUj(Z''- 6*^) into itself. So {TZ'l;{a-9j)\n e N} 
is a finite set. Therefore there exists some n G N, and q > such that TZ'^{a — Oj) = TZ^''{a — 6j)- Thus 
TZ^{a — Oj) is periodic. 

From (i) we have that —TZ^{a — 9j) is cyclic for {Td)dev- So ~TZ'^{a — 9j) is in the attractor X{A^ ,T>). 

(iii) From (i) and (ii) it is clear that the periodic points for TZc lie in {U^ — C) fl {—X{A'",V)). For the 
other inclusion take ti G {U^ - C) n {-X{A^,V))- Then using the formula (|6.3p , there exist di,d2 - ■ ■ & T> 
such that 

-ti = (^^)-idi + (A^)-2d2 + ... 
Let -t„ := (A'^)-M„ + (A'^)-2d„+i + . . . . We have 

(6.11) A^(t„) + d„ -t„+i, (nGN). 

Since ti G Z'' — C, equation (16. lip implies that t2 is in U^ — C and ti — TZcit-i)- By induction t„+i is in U^ — C 
and 7^c(^,^+l) = t« for all n G N. But we have also t„ G -X(A'^,I?). And, since (Z'' - C) n (-X(yl'^, X>)) 
is finite, there exist n, m > 1 such that t„ = tn+m- This implies that 7^™(t„) = TZ^{tn+m) = t„. Since 
7?.p^^(t„) = ti, it follows that ti is periodic for TZc- D 

Theorem 6.15. Let C := {6'o, • . • , ^'p-i} &e a simple cycle. 

(i) For eac/i fc G Z'' bkc? eac/i j G Zp t/iere is a unique uo{k^i) — lo^LjOi • • • G 51 such that for all x G , 
(6.12) 

Moreover there exists a cycle C ^ [C — U^) n X(A^,I?) smc/i t/iat u!{k, j) G ilc'- 

(ii) r/ie infinite word uj{k, j) can be constructed as follows: 

(6.13) ujn^Tll.{k-ej)-A^TZl+\k^ej), (n G N). 

(iii) Suppose C is a cycle in {C — Z'^) O X{A^ ,T>), and uj G ilc" - Then there is a unique {k{uj) , j {uj)) G 
Z'^ X Zp SMc/i that for all x G R"^, 

(6.14) (g2.,.^g2..r.„x^g2..r.,r.„.^ ^ ^ ) ^ ^^(^ ^ _ e,M,j('^))- 

(iv) {k{uj) , j {ll})) can be constructed as follows: if lo Cz flc has the form loq . . - LUnp-imo . . .m,g_i, then 
the fixed point rjo of t„i^_i - - - Tm„ belongs to Oj[^) — Z'' for some unique j{Lu) G Zp. And 

(6.15) k{u;) = c^o + • • • + (A^)"P-^c^„p-i + - (A^r^vo- 
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(v) Let 

The maps 
and 



nc := |J{^C" I C" cycle in (C - ) n X{A^,V)}. 
ec : Z"* X Zp ^ fie, ^cikj) = j), 



are inverse to each other. 

Proof. Let (k.j) e Z'' x Zp and let uj(k,j) be defined as in (ii). We prove that the relation (|6.12p is satisfied. 
We have 

A^TZcik - Oj) + Wo = fc - Oj^A^TZlik - 9^) + wi = Ucik - Oj), . . . 

Therefore 

7^c(fc - 6j) = {A^)-\k - e^) - {A^r^u:Q,nl{k - dj) = {A^)-\k - 9,) - {A^y^cjo - {A^r^u;,, . . . 
By induction 

So 



T^cik - e,) = (A^)-"(fc - e,) ~ {A^r^LUo [A'^r'u;,,^ 



r^„_i . . . r^ax = [A^r^x + {A^r^uo + ■■■ + {A^)-'lo,,^i = + iA^)-^{k - 6,) - ni{k - 6,). 

But 7^^(fc -Bj) <E - 0j+n so 

T^„_i ■ • ■ T^oX = (^^)""(a; + k- Oj) + Oj+n- 

Therefore the relation (|6.12[) is satisfied. 

Next we prove the uniqueness of uj. Suppose lo' also satisfies (|6.12p . Then t^'^x = r^gX so {A'^)^^uJq = 
{A^'j^^ujQ which imphcs that uj'q — ujq e Since I? is a complete set of representatives for 11^ j A^U^ , it 

follows that cjg = Wo- By induction = a;„ so w' = w. 

To see that uj{k,j) is in some ujc' for a cycle C in (C — Z'') n X{A'^,'D) we use Proposition [OH There 
exists an I such that 'R}(j{k — 6j) is periodic for TZc, so — 7?.^(fc — 9j) is a cycle point for {Tii)dev- Let C its 
corresponding cycle. By Proposition 16. 141 C is contained in (C — Z'^) n X{A^ .V). By the construction of 
Ljj{k,i) given in (ii), and by Proposition 16. 14f i) . we see that uj{k,p) G Vlc' ■ 

Now let w be of the form given in (iii). And let {k,j) :— {k(uj), j{u!)) be as in (iv). Since rjo £ 9j — Z''' it 
follows that {A'^y^Prjo is also in 9j — Z'' so fc = k{u}) is indeed an integer. 

We check that if for (fc, j) we construct f = i/qUi... = uj{k,j) as in (ii), then v = uj. This will prove also 
(v). From (I6.15P we have with to = —{k — 9j), 

h ■■= T^oito) - T^o(-(fc - Sj)) = -UJi {A^rP-^Unp-l + {A^r^-h^o e 9, + i - Z\ 

since 770 e 9j - Z''- implies that {A^Y'p-'^tjo e 6'j-(„p-i) -Z'^ = 9j+i - Z'^. This shows that 7^c(-^o) = -^i, 
and V{) — L0{). By induction we obtain vi = wi, . . . , Vnp-i — ^np-i and that —rjQ = TV^^—to). And since the 
cycle C of 7/0 is given by mo . . . rriq^i, it follows by Proposition 16. 14^ 1) that j/ = wo . . . ujnp~imo . . . viq-i. 
For the uniqueness of (fc, j), suppose {k,j) and {k' ,j') satsify (|6.14[) . Then 

(A^)-"(fc - 9,) + 9n+j = {A^)-'\k' - 9,,) + 9,,+j,, {n e N). 

Taking the limit as Tip 00 we obtain 9j — 9ji E Z'^. Since the cycle is simple, j = j' . Therefore 

{A^y{k - 9j) = {A^y\k' - 9,), {n e N). 

But this means that i{k — 9j) = i{k' — 9j), and by Proposition 14. lU is injective so k — k' . □ 

We summarize our results in the following corollary. 

Corollary 6.16. Suppose (A,!)) satisfy the tiling condition fDefinition \6.6]) . Let C — {9o, ■ . ■ ,9p-i} be a 
simple cycle. Let 

nc := \J{nc' I C cycle in [C - Z'^) n X{A^ , V)}. 
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(i) The maps 

d:X{A^,V) X nc^SAiC), d{x,uj) = (e2--, g^^-^o-, e'^'-i'-o", . . . ), 
ic:R''xZp^ SAiC), icix,j) = (e2"«-4")""^+''"+^-))„eN 

are bijections. 

(ii) 

where k{u!),j{Lu) are defined in Theorem 1 6. 1 5Y iv } . 

d-\ic{x,j)) - {y,io{k,j)), (x e R^j e Zp), 

where y £ X{A'^ ,T)), k £ 1/^ are uniquely defined by x + 9j = y + k, and uj{k,j) is defined in 
Theorem \6.15Y ii ) . 
(iii) The following diagram is commutative: 

X{A^,V) X Vtc — °— ^ Sa{C) ( R'' X Zp 



X{A^,V) X Vie — °— ^ Sa{C) ( M'^ X Zp 

Corollary 6.17. M^zi/i the notations in Theorem \6. 151 we have 

(6.16) j{uJiUJ2 ■ ■ ■) ^ jiujQUJi ■ ■ .) + I, {uJauji--- eftc). 

(6.17) (A^)""(a; - 6'j(^^o'^i--- ) + ^(^owi ■■.))= t^„_i ■ ■ . t^o^ ~ Sj(<^oi^i... )+« + fc(cj„tJ„+i . . . ) 
/or aZ/ X e X{A^ ,'D),LUoUJi ■ ■■ £ Clc- 

Proof. We apply the commutative diagram in Corollary 16 . 1 61 to p": 

i'^Hp^{x,uj) = «c5(7"(^„_i • • • T^aX, w„a;„+i . . . ) = (Ti^„_i . • ■ T^(,2;-^'i(i^„i^„+i... )+fc(WnWn+l • • ■ ), j(a;„Wn+i ...)). 

"a!^^c^^(2^'^) = ~ ^jX'^oc^i...) +fc(wowi.-.):j(^oWi---)) = 

((A^)-"(a; - 0j(-oc.i...) + fcl'^owi ■ • ■ )), jX^o^^i • . . ) + n). 
Since the two quantities are equal to each other (by the commutative diagram in Corollarv l6.16p . the relations 
follow. □ 

With the aid of our cycles and associated encoding/decoding mappings we are now able to state our 
main result regarding super representations. Notice that the introduction of cycles yields the following 
improvement of Theorem 14.91 in section 4 above. 



Corollary 6.18. Suppose {A,T>) satisfies the tiling condition, and let C be a simple cycle of length p. On 
X{A^,T>) X (ic define the measure fi by 

I fdfi^ I V fix,uj)dx. 

Jx(AT.V)xflc Jx(AT,V) ~ 

Define the operators fk, k £ Z'^ and U on L'^{X{A'^ ,T>) x ^ci^) by 

fkfix^uj) = e2-*-/(a;,c.), (x £ XiA^,V),cj £ ilc), 



Uf^^/\d^\fop-\ 

Then {fk, U} define a unitary representation ofGA andW : L'^i'R'^xZp) L'^ {X {A^ ,V) , ji) , Wf = fodoi^^ 
is an isomorphism that intertwines this representation with the one in Definition \4-.6\ 

When {A, T>) satisfy the tiling condition we can say a bit more about the possible extra cycles in (C — 

Z'^)r^x{A^,v)■. 

Proposition 6.19. Suppose {A,T>) satisfies the tiling condition. Assume that there is a cycle point 0q £ 
X{A^,T>) such that Oq — k £ X{A^,T>) for some k £ Z'^ , fc ^ 0. Then the entire cycle of Oq is on the 
boundary of X{A'^ ,T>). 
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Proof. Let denote the interior of V). We will prove first that if a point x is in X{A'^, V) n 

{X{A^ ,1)) + k) with k e Z'', fc ^ 0, then x is on the boundary of X{A^,'D). Suppose not, then 
X e X{A^,V)°. By jLW96c| we know that the closure of X{A'^,V)° is X{A^,V). This implies that 
the neighborhood X{A'^,V)° of x must intersect the set X{A^,V)° + k. But since ^(A^,!?) tiles R'^ by 
Z"*, this implies that the interiors of X{A'^ ,1)) and X{A^,'D) + k cannot intersect (the intersection would 
have positive Lebesgue measure). So x must be on the boundary of X{A'^ ,1)). 

Now consider 9q and let C = {Oq, . . . , 9p-i} be the cycle of 9q and let Iq, . . . Ip-i, the corresponding digits. 
We have ti^_-^ . . . n.Oj = for all g {0, . . . ,p — 1}. Suppose one of the points Oj of the cycle C is in 
X{A^,V)°. Since 

X{A^,V)= y r,(X(A^,I?)), 

dev 

we obtain that Td(X(A-'", X>)°) C X{A^ ,1))° {tcl is a homeomorphism) . So if 9j is an interior point for 
X{A^ ^V), then Oq = ti^_-^ . . - Ti-Oj is also in the interior of X{A^ ,1)). This contradiction implies that C is 
contained in the boundary of X[A^ , V). □ 

To help the reader appreciate our encoding results we present some examples which at the same time 
stress tiles versus spectrum. Since the technical points are illustrated already for the real line we begin with 
dimension one, and then turn to the plane M^. 

Example 6.20. Let us take d — 1, A — 2 and V — {0, 1}. The maps are tqx — x/2, tix = {x + l)/2. 
Consider the simple cycle C := {9o} = {0}. It corresponds to 0. 

The attractor X{A'^,V) is [0, 1]. The intersection (C - Z) n X{A^,V) = (-Z) n [0, 1] = {0, 1}, so it 
consists of the cycles C": {0} and {!}, which correspond to and 1 respectively. Therefore 17o = f^o U ili, 
i.e., the words that end in an infinite repetition of or an infinite repetition of 1. 

We have the map 

1(7^0 : [0, 1) X (Qq U Qi) R, uj) ^ x + k{uj), 

(we used [0, 1) here instead of [0, 1] to have that the map iQ^d is a true bijection, not just up to measure 0), 
and with formula (|6.15p : 

k{uJo . . . ujnO) = tJo + 2 • cji H h 2"uj,i, K^o ■ ■ ■ ^ni) = Wo + 2 • t^i H h 2"tj„ - 2"+^ 

Applying the commutative diagram in Corollarv l6.16l to p", we have iQ^dp"'{x,uj) ~ ■^iQ^d{x,uj), which 
implies that 

■^{X + k{uJoLUi ...))= T^„_i . . - T^gX + k{u!nUJn+l ■■■), {x G [0, 1) , U! £ Hq U fti). 

Example 6.21. Let d = 1, A = 2 and I? = {0, 1}. Consider the simple cycle C := {6*0, 6*1} = {1/3,2/3}. 
It corresponds to 10 (because ti(1/3) = 2/3, ro(2/3) = 1/3. The attractor is X{A^,V) [0,1]. Then 
(C - Z) n [0, 1] = C, so Clc = ^c, i.e., the words that end in 10. 
We have that the map 

l^^d : [0,1) X 17c X ^2, ic^Hx^u;) = (x - 6lj(„) + fc(w), j(w)) 

j{uJo ■ ■ ■UJ2n~llQ.) = 0, j(^o • ■ .a;2n-i01) = 1, 
k{LUo . . .C^2„-ll0) = C^o + 2 • C^l + • • • + 2^"-^LU2n-l + 2^" • 1, 

fc(wo • • • W2„-i01) = c^o + 2 • c^i + • • • + 22«-icj2„-i + ^ - 22" • I 
The map ftc 3 uj i-^ {k{uj), j{Lu)) e Z x Z2 is a bijection. 



is a bijection, and 
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As an example, let us show how to compute the uj G fie associated to (fc, j) — (15, 0). Take k — Oj = 15— -I- 
We want a fci g Z and G 2? such that 15 — i = 2(fci — §) + wq. We have 





1 


= 2(8- 








1 








15- 


3 


5) 


+ 


1 - 


3 


= 2(1- 


5) 


+ 




2 


= 2(4- 
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= 2(0- 






8- 


3 


5' 


+ 


1 - 


3 




+ 1 




1 


= 2(2- 
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= 2(0- 






4- 


3 


1' 


+ 1 


0- 


3 




+ 1 




2 


= 2(1- 








2 


= 2(0- 






2- 




^) 


+ 


0- 






+ 




3 








3 







Thus cj(15,0) = 00100110. 

Remark 6.22. Our next example is in the plane, but it illustrates a more general picture in M.'^ for any d. 
Start with a given pair {A, V) with the matrix A assumed expansive, and 2? a chosen complete digit set, i.e., 
in bijective correspondence with the points in Z'^/A'^Z'^. So in particular, = | det A|. In general it is not 
true that the same set I? is a digit set for A, i.e., that it is a bijective image of Z'^/AZ'^. Here for d = 2, we 
give an explicit geometric representation of a pair {A, V) for which the same I? is a digit set for both the 
radix representation with A and with the transposed matrix A^ . Hence we get two attractors X{A^ ^T)) and 
X(A, P). Both will be referred to as Cloud Nine, a left-handed version, and a right-handed version. That 
is because there are nine integer points, i.e., the intersections with consists of nine points, and it is the 
same set for the two fractals. For each, there are three one-cycles, and one six-cycle. While the six-cycles 
(for A and for A^) are the same as sets, we will see that they are traveled differently under the actions 
discussed in our encodings from sections 5 and 6 from above; the difference being essentially a reversal of 
orientation. Hence our encoding with infinite words in letters from T) will also be different for the two cases, 
and the details are worked out below. Recall the attractor X{A^ ^V) is the set of "fractions" for our digital 
representation of points in . The attractor X{A^ ^V) is also an affine Iterated Function System (IFS) 
based on [A^V). Thus the Cloud Nine examples further illustrate the intricate part played by the cycles in 
Z^ for the initial {A, r')-IFS. In each case, using these cycles we are able to write down formulas for the two 
maps which do the encoding as well as the decoding in our positional P-representation. 

For use of matrices in radix representation, the distinction between the radix matrix A and its transpose 
A^ is important. First the two matrices sit on separate sides in a Fourier duality; and secondly, even if the 
chosen set of digits V is the same, the two attractors may be different. In fact, in general the same V may 
not work for both A and . There is not a natural connection between the two quotients Z'^ jAZ'^ and 
jji- j X^jd. ^ J ^ ^T^^ other one for the transposed. 

But in the particular 2D example below. Example I6.23[ called Cloud-Nine, one may check by hand 
that, for this matrix A, with det A = 5, each of the two quotients Z? j Al? and Z? jAJ^l? are in bijective 
correspondence with the same subset V m.Z? . (See details!) 

As a result it makes sense to analyze the two different Hutchinson attractors X(yl^,I?) and X(A, P); 
both compact and with non-empty interior. The first one, in a different context was studied earlier in [BJ99] 
and |Jor03| . but both are interesting. Note that while two references |BJ99j and |Jor03| use these examples 
the questions addressed in these papers are completely different. 

As we see, there is an intriguing connection between cycles, solenoids, and encodings for the two. 



Example 6.23. Take d = 2, A = ( ^ ) ( 1 ) ' '^'^'^ ^ 








±3 












7 




±2 



We consider the trivial cycle C := { 



}. We want to compute (C -Z'^)r\X{A\V), i.e., Z'^^X{A\V) 





^ 

First, we need to locate the attractor X{A^,'D). For this we use the proof of Proposition I6.14r ii). and 
conclude that if R := ^^^^ i-h(a^ ' )^-i'|" ^^'^^^ then the ball B{0, R) is invariant under all the maps Td which 
implies that X{A'^,'D) is contained in this ball. 

Since ||(A^)-i|| = we conclude that R = = 2.427. . .. There are 21 points in Z'' n 5(0, i?) and 

we can check how TZc acts on each of them. 
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If we want 
(6.18) 

Thus, given 



A' 



di 
d2 



with 

{x - d,) - 2{y - d^) 



e 1? and [ 

"2 

2(a;-di) + (y 



G I?, then we must have 



to find Tlf 



X 




a 


and 


■ rfi " 


. y . 




b 


. ^2 . 



e T), we first look for an element in T) with 



di — 2^2 ^ X — 2y mod 5 and 2di + d2 = 'ix + y mod 5. Then we compute a, b as in (|6.18p . 

It is interesting to note here also that if di — 2^2 = x ~ 2y mod 5 then the other equivalence mod 5 is 
satisfied too. This is because 2di + d2 = 2{di — 2^2) mod 5 etc. We can use the following table 





di 













3 






-3 




















d2 
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-2 




(di- 


2d2) mod 5 
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2 


1 


4 



We apply these ideas to the points in 7? O B{0, R): 
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+ 
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■ 1 ■ 


+ 








so TZc 

so TZc 



-2 



" ■ 









" 1 ■ 










so TZc 



and — 
and 



is a cycle that corresponds to 
is a cycle that corresponds to 



and 



is a cycle that corresponds to 
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■ -1 ■ 
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■ 3 ■ 




■ -1 ■ 
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1 ■ 
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■ ■ 
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-1 
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-1 
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' 1 ■ 




■ -3 ■ 




■ 1 
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1 



= AT 



1 ■ 




■ 




■ -1 " 
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-2 
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1— ' 



= A^ 



1 

-1 

1 



= A^ 



■ 




CO 
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-1 


+ 
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-1 



-3 




So 






TZc 


-1 


Uc 


1 







Tic 
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TZc 


-1 







1 




-1 




-1 




-1 




1 




1 




1 



and 





1 




-1 









-1 




1 
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1 




1 




-1 




-1 



is a cycle that corresponds to the word 



CO 
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-3 









-3 
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-2 








Similar computations show that 
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TZc: 
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TZc 
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TZc 
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1 



TZc 





-2 


TZc 


-1 


TZc 












-2 




-2 



-2 


7^c 







-2 







1 




-1 
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-1 




-1 



TZc 



TZc 



Thus we have three cycles of length one, and one cycle of length 6. 
Consider now the pair [^A^ , T>) , that is replace the matrix A by A^ , and T> is also a complete set of 
representatives for 'Z,'^/{A'^)'^Z'^. As above we get the same cycles: 



is a one-cycle that corresponds to 
is a one-cycle that corresponds to 
is a one-cycle that corresponds to 



We also obtain the six-cycle 
to a different periodic word: 



1 

-1 



-1 
-1 





-1 



-1 
1 



This cycle corresponds 
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-3 









-3 







-2 
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This means that the 6-cycle is traveled by two diferent paths according to the matrix A or . 

6.2. Tiling and spectra in some examples. One of the uses of encoding is apphcations to tiling questions 
in M!^. The simplest tiles X in M.'^ are measurable subsets which make a tiling of M/^ by translations with 
vectors from some lattice, say F (i.e, a rank d subgroup) . Since we work in the measurable category we allow 
different translates in the tiling X + 7, for 7 G F to overlap on sets of measure zero. 

One might think that when X is given, then the presence of a suitable lattice F making X into a translation 
tile for M.'^ could be decided by visual inspection, at least in the case of d = 2. After all, when a pair {A, V) 
is given, then there are fast Mathematica programs which produce excellent plots of the attractor sets 
X = X{A'^ , V), black on white; see for example [B J99| . But except for isolated cases, it turns out that when 
the planar sets X are represented in black on white, then there will typically be many white spots, or gaps, 
disconnecting X in complicated ways. If some lattice F will make X into a translation tile, then the white 
areas must be filled in by black under translations X + 7, 7 € F. An inspection of [BJ99| reveals that this 
is not easy to discern by visual inspection. Hence, instead we resort below to spectral theoretic tools for 
locating the lattices which do the job. 

A more complicated form of tilings (still with a single base tile) refer to the case when the set F of 
translation vectors is a set which is not a lattice, e.g., translation sets of quasiperiodic tilings. But for our 
present considerations lattice tilings will suffice. 

The sets X which will interest us are the attractors X = X{A'^,'D) from afhne IFSs as described in 
CoroUarv 16.161 It is known that every such X is compact with non-empty interior, and so in particular it 
has positive d-dimensional Lebesgue measure. 

Hence it is of interest to ask for a spectral analysis of the Hilbert space L'^{X), referring to d-dimensional 
Lebesgue measure. In fact, using Pontryagin duality for abelian groups, one can check that X tiles with 
a lattice if and only if the dual lattice makes an orthogonal basis of complex exponentials in L'^{X). The 
result is often refered to as Fuglede's theorem. For background, see the references |Fug74| and |Rud62| . 

To understand the correspondence between choice of translation lattice on the one hand and spectrum on 
the other we need: 

Definition 6.24. Let X C R'^ be measurable with < niX) < 00 where fi denotes the d-dimensional 
Lebesgue measure. For ^ e M'' set e^{x) = e^'^*^'^, where ^ • x :=< ^, x >= ^iXi + ■ ■ ■ £,dXd and x = 
{xi, - ■ ■ , Xd) & M''- If F C M"* is a discrete subgroup (in this case a rank d lattice) set 

Ex{r) :={ee|x : C e F} 

where \x denotes restriction to the set X. 
If F is a lattice we set 

F° := {A e R'' I A • ^ e Z for all ^ € F} 

called the dual lattice. 

If A C K'' is a discrete subset we say that it is a spectrum for X or that the pair (A", A) is a spectral pair iff 
Ex {A) is an orthogonal basis in the Hilbert space L'^{X) = L'^{X,pi). 

Lemma 6.25. (Fuglede |Fug74| ; Let < ^l{X) < 00 and A a rank-d lattice. The following conditions are 
equivalent: 

i) Ex {A.) is an orthonormal basis in L'^{X); 

ii) X tiles by the dual lattice A°. 

Remark 6.26. We can draw the following stronger conclusion: When X is given, there are no other tiling 
lattices for X than those which arise as in (ii) by spectral duality. Proof. Every lattice A satisfies A°° — A, 
i.e., the double dual yields back the initial lattice. To see this, use the following general observations which 
also serve to make explicit the standard lattice operations which we will be using in the proof. 

Referring to the IFS of Dcfinition l6. ll we note the following formula for the computation of the L'^{X{A'^ , T>))- 
inncr products. Set X = X{A'^,'D) and Xx{£,) ■— fx e^{x)dx. Then 

00 

Xx{i) = \{mv{A~^0 

n=l 
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where 



1 



|detA| 



d<£V 

Recall that | det A\ =number of elements in V. 

Some remarks about lattices in K'' are in order: every lattice is by definition a rank-d subgroup of R"* 
and it can be shown that it has the form F = M'E'^, where M is an invertiblc dx d matrix and where points 
in Z'' are represented by column vectors. We will write Tm to emphasize the matrix M that completely 
determines the lattice. The next lemma is elementary: 

Lemma 6.27. (i) Tm C Z''- if and only if M e Md{Z)- 

(ii) IfT = Vm then V° — r^^^T^-i. In other words ifV is given by M then its dual is given by {M'^)^^ . 

(iii) r°° = r 

We will use names from |BJ99) for the fractals X = X{A'^,'D) in R^. These names refer both to their 
geometric appearance as planar sets X, as well as to a counting of Z^-cycles, i.e., the number of points in 
{~X) n Z^. See |BJ99| (end of subsection 9.3) for details. For example, Cloud-Nine has three one-cycles and 
one six-cycle in Z^. 

What follows is a family of examples in 2D. In each case, we are asking the following questions: How 
much flexibility is there in selecting digits when the base for the 2D number system is fixed? In our case, 
we are using the positional radix representation for vectors, and thus the base for our number system is a 
chosen matrix A. For several of the examples below, we fix a particular A, and then we vary our choices of 
"digit" sets I? in Z^. The points in V will serve as "digits" in a positional representation. 

We are motivated by Knuth's algorithmic approach mentioned in the Introduction: What are the "inte- 
gers" and what are the "fractions" in a number system specified by a particular pair {A,!))? What is the 
encoding, and what is the decoding? When the matrix A is fixed, how do changes in 1) reflect themselves in 
the answer to the questions? 

The examples below are sketched with Mathematica programming in |BJ99j . and the names we use for 
the fractals X are consistent with |BJ99| . i.e., the the respective names of the sets X, Cloud-Nine etc. The 
examples when A is the same but V changes are referred to by the name Cloud, followed by a number. The 
number indicates the cardinality of {—X) n Z^. 

However the questions addressed here are different from those of [B J99] . 

It is of interest to understand how much flexibility there is in selecting digits when the base for the number 
system is fixed. In our case, the base for our vector number system is the matrix A, and so we vary the 
choices for the companion set P. But when A is given, the choice of D is always restricted by demanding a 
bijection V ^ 1?IJ(^1?. 

' 1 -2 



Cloud Three. A 



V 



2 1 

tile K'^. Cloud Three only has one-cycles on Z 
wavelet. It has measure — 2. 

1 -2 

2 1 

lattice Z x 2Z. So X is not a Haar wavelet 

1 -2 

2 1 


















±1 





i-x) 




±2 

X) n z' 



Cloud Five. A 



Cloud Nine. A = 



V 



V = 








±3 







7 





7 



±1 



It has measure = 2 








±3 















7 


±2 



the lattice Z x 2Z. Cloud Nine X has three one-cycles and one six-cycle 
Haar wavelet. It has measure — 2. 



Here the lattice 2Z x Z makes X 
= Ci. Moreover X is not a Haar 

Cloud Five is a lattice tile with 

Cloud Nine is a lattice tile with 
So Cloud Nine is not a 



Twin Dragon. A = 



1 



-1 



V = 



The Twin Dragon is a lattice tile with 



lattice Z'^ = F. So the Twin Dragon is a Haar wavelet. It has measure — 1. 
We are using Lemma 16.251 in identifying lattices which make the various Cloud examples X tile . For 
this purpose we must identify our cycles relative to so called Hadamard systems as defined in |DJ07cj . A 
Hadamard system consists of a matrix A and two sets T> and £ as dual digits, = = \ dctA\. By 
"dual" we mean that the matrix formed from the exponentials as 

(6-19) — ===(e ^ ' jdeVMC- 

Vidct A| 
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is a unitary | det A| x | det A| matrix. 

Let N be the absolute value of the determinant, and let Zjv be the cyclic group of order N. Then the 
matrix Un for the Fourier transfrom on Z^r is an example of a Hadamard matrix as in (|6.19p : specifically 
the j, k entry in Un is "^C"'''; ^ ^ '^N, where C = Cjv is a fixed principal A^'th root of 1. 

Proof. To find the lattices that give tiles for these examples, we use Lemma [6.251 and find the dual lattices 
that give orthogonal bases of exponentials. For this we use the techniques introduced in [D J06b[ [D J07c| . 

/ 1 -2 \ 

First let us look at the matrix A — [ for the cloud examples. We want to find what is the 



lattice that makes X{A'^,'D) = \dj e V} tile . For this we need a set C such that 

1 u2iTt(A^) ^"i '^^gp is a unitary matrix, i.e., (A^, V, C) is a Hadamard triple. It is enough to take 

£ a complete set of representatives for Z^/AZ^. We will take C to be 



C 








±3 
















±2 



With this choice of C the reader may check that for each of the Cloud-examples listed above, the corre- 
sponding Hadamard matrix from (|6.19p turns out, up to permutation, to simply agree with the matrix U5 
for the Fourier transform on Z5. 

According to |DJ07cj we have to see if there are any proper invariant subspaces for A. But those would 
give rise to real eigenvalues of A, and this is not the case. Thus, by [DJ07c| . the spectrum of X{A'^,'D), is 
determined only by the "m-p-cycles" . These are the cycles C — {xq, ■ ■ ■ ,Xp-i} for the "dual" IFS a-i{x) = 
A~^{x + I), I £ C with the property that \mD{xi)\ = 1 for all « S {0, . . . ,p — 1}. 

Then, by [DJ06bl lDJ07c| the spectrum of X{A'^,D) is the smallest set A that contains — C for all the 
TO-p-cycles, and such that AA + £ C A. 

Cloud Three: We have 

mv(x, y) = -(l + e^'^'y + e^^^^y + g^'^^^y _^ g-27^^2^;^ _ 
5 

If we want |mx)(x, y) \ ~ \ then we must have that all the terms in the sum are 1 so y e Z, and x is arbitrary. 
We are looking for mp-cycles, so |mp((T/(a;, y))| — 1 for some ^ G £ so ai(x,y) must have the second 



component in Z. The inverse of A is A ^ — i ( ^ ) . Thus ■^{2{x + l,) + (y + ly)) G Z. This implies 

that X G 5Z. 

We claim that A = ^Z x Z. For this, note first that A{^Z x Z) +£ C ^Z x Z. By the previous computation. 



2^^N^. uiiio, iii^^^ iiiou ^i^2'"'^'"; ''-^2 
iZ X Z contains the negative of all the mp-cycles. Then, take wq '■— (|, k') G ^Z x Z. Then for I G £, 

ai{-^,-k') - i^i-k + 21, + 4k' - 4ly), + 21, ~k' + ly)). 

Note that there is a unique Iq E C such that wi :— ~o-ig{—^, ~k') G ^Z x Z. As in the proof of Proposition 
16.141 there is a sequence IqJi - ■ ■ G £ such that if Wn+i = —Ti„{~Wn), n G N, then w„ G ^Z x Z and, for 
some m, —Wm is a cycle point for {(ti)i^c- Note that since Wn G ^Z x Z, —Wm is a point in a mp-cycle. 

Since Wm = —Tim-ii'^m-i), we have that Wm-i — Awm + Im-i G A{—C) + £, where C is the m-p-cycle 
of —Wm- By induction we obtain that wq must be in A. Thus A = ^Z x Z is the spectrum. Taking the dual 
we obtain that X{A^ ,V) tiles by 2Z x Z. 

Cloud Five: We have 

mnix, y) = -{l + e^^^*^^ + g-^'^^^x ^ ^2^^x ^ ^-2^tx\ _ 
5 

Therefore \mT>{x,y)\ = 1 iff a; G Z. For mp-cycles we must have that the first component <7i{x,y) must be 
in Z for some / G £. This implies that y G ^Z. 

We claim that A = ^Z x Z. The proof works just as for the Cloud Three example so we will leave it to 
the reader. Thus the dual lattice is 2Z x Z and X{A'^,'D). 

Cloud Nine: We have 

mv{x,y) = - (1 + e^^^*^^ + e-2^*3:r _^ ^2rrz2y _^ ^-2^^2y^ 

5 
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So \mTi{x,y)\ = 1 iff .T e iZ and y € ^Z. For mp-cycles we must tiave {x,y) = (|, ^) with k,k' & Z and 
also |m25(crz(|, -y))! = 1 for some ^ G £. This implies that ^(2(| + dx) - + dy)) G 5Z so fc must be 
divisible by 3. Thus the m-p-cycles are contained in Z x ^Z. 

Just as in the previous examples we get that A = Z x iZ so the tiling lattice is Z x 2Z. 

Twin Dragon: For A = ^ l ^ ' ^^^^^ ^'^ proper invariant subspaces so the analysis of the 

TO-p-cycles will suffice. We can take £ := 

mvix,y) = ^{l + e'-^^). 

Therefore \mxi{x, y)\ = 1 iS x E Z. A^^ ~ 5 ^ | 1^ ) ' want the first component of ai{x, y) to be in 

Z so — y) e Z, therefore y G Z. 

As in the previous examples we can check that A = Z^ so the dual lattice is Z^. □ 

The method used in our analysis of the examples may be formalized as follows. Stated in general terms it 
applies to a large class of IFSs which carry a fairly minimal amount of intrinsic duality. For the convenience 
of the reader, we begin with two definitions. This is of interest as there are few results in the literature 
which produce formulas for lattices which turn particular attractors X into tiles under the corresponding 
translations. For details we refer to |DJ07c( fDJOSbj . 

Definition 6.28. (i) A Hadamard triple in is a system (A,!?, £) where A is expansive in Aid{Z), 
are in Z"^, and C is such that the matrix (j6.19p is unitary. 

(ii) For a Hadamard triple {A^V.C) the cycles C, that correspond to the IFS {cri)ieC: for which the 
absolute value of is 1 are called extreme relative to m-p, or mxj-cycles. 

We are now ready to state our general tiling result. 

Corollary 6.29. Let (A^DjC) be a Hadamard triple in W^, and such that D is a complete set of represen- 
tatives for Z^/A'^Z'^, and let X = X{A'^,'D) be the corresponding Hutchinson attractor. Suppose A has no 
proper invariant subspaces. Let A be the smallest lattice in M.'^ containing all the sets — C where C runs over 
the m-p-extreme cycles, and which is invariant under the affine mappings x 1— > Ax + I, for I G C Then the 
dual lattice F = A° makes X tile with T translations. 

Acknowledgements. The second named author had helpful discussions with Prof Sergei Silvestrov, University 
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